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ABSTRACT 

We demonstrate that the Biermann battery mechanism for the creation of large 
O \ scale magnetic fields can arise in a simple model protogalaxy. Analytic calculations and 

q \ numerical simulations follow explicitly the generation of vorticity (and hence magnetic 

field) at the outward-moving shock that develops as the protogalactic perturbation col- 
lapses. Shear angular momentum then distorts this field into a dipole-like configuration. 
The magnitude of the field created in the fully formed disk galaxy is estimated to be 
1CP 17 Gauss, approximately what is needed as a seed for the galactic dynamo. 

o 
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i-^ 1 1. Introduction 

The origin of galactic magnetic fields has proved to be one of the most challenging and stubborn 
problems in modern astrophysics (Rees 1987; Kronberg 1994 and references therein). It is generally 
assumed that galactic fields are generated and maintained by the dynamo action of a differentially 
rotating disk galaxy. However, a dynamo can only amplify an existing field and so the question of 
^ ■ galactic magnetic fields splits naturally into two parts: creation of the field required to seed the 

dynamo, and the nature of the dynamo itself (e.g., Zel'dovich, Ruzmaiken, & Sokoloff 1983). 

An early attempt to explain the origin of seed fields is due to Harrison (1970, 1973) who 
showed that magnetic fields are created during the radiation era if significant vorticity exists at that 
epoch. However, primordial vorticity in an expanding universe decays with time (in contrast with 
the irrotational density perturbations presumably responsible for structure formation). Indeed, 
the absence of significant vorticity prior to galaxy formation together with the observation that 
vorticity is generic to galactic disks may provide an important clue as to the origin of galactic 
magnetic fields. Angular momentum in galaxies is thought to arise from tidal torques among 
neighboring protogalaxies (Hoyle 1949, Peebles 1969, and White 1984). However, gravitational 
forces alone do not produce vorticity and therefore its appearance must be due to 'gasdynamical' 
processes such as those that occur at oblique shocks. These same processes also produce magnetic 
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fields by creating a so-called Biermann battery (Biermann 1950) which drives electric currents in 
the plasma. Since oblique shocks are inevitable in collapsing gas clouds, the early stages of structure 
formation provide a natural site for the production of seed fields (Pudritz & Silk 1989, Kulsrud et 
al. 1997). 

There have, of course, been other attempts to explain the origin of seed fields. One possibility 
is that first fields were created in stars and subsequently expelled into the interstellar medium 
(Bisnovatyi-Kogan, Ruzmaiken, & Syunyaev 1973). Alternatively, seed fields may have been created 
in the very early Universe through such exotic phenomena as quantum field creation during inflation 
(Turner and Widrow 1988, Ratra 1992), phase transitions (Quashnock, Loeb, & Spergel 1988, 
Vachaspati 1991, Field & Carroll 1998), and topological defects (Sicotte 1997). By comparison, 
the protogalactic battery has a certain simplicity and elegance since the necessary ingredients are 
generic to models of galaxy formation. 

Most discussions of galaxy formation ignore the influence of magnetic fields (see, however 
Wasserman 1978 and Kim, Olinto, & Rosner 1996). To be sure, galactic magnetic fields play an 
important role in a number of astrophysical processes such as star formation, cosmic ray confinement 
and gasdynamics. But while the energy density in galactic magnetic fields is comparable to that in 
cosmic rays and in the turbulent motion of the interstellar medium, it is considerably less than the 
energy density associated with the global dynamics of a galaxy. This suggests that magnetic fields 
play a secondary role in the formation and evolution of galaxies. Nevertheless, models of galaxy 
formation should be able to explain their origin. 

In this work, we investigate the generation and early evolution of vorticity and magnetic fields 
in the context of a detailed, albeit highly idealized, model protogalaxy. Recently, Kulsrud et 
al. (1997) have attempted to follow the creation of protogalactic magnetic fields in a cosmological 
hydrodynamic simulation of a cold dark matter universe. They find that seed fields can be produced 
on a variety of cosmologically interesting scales. Our work complements and enhances theirs by 
considering a simpler system where both analytic and numerical techniques can be employed. In 
so doing, we are able to understand these earlier results and some of their limitations in terms of 
relatively simple physics and numerics. Our work also makes contact with various semi-analytic 
models of disk galaxy formation. 

An outline of our scenario is as follows: 

• We consider an isolated, nearly spherical, density perturbation in an otherwise Einstein-de 
Sitter universe. The cosmic fluid consists of both collisional gas and collisionless dark matter. 
We assume that the scales of interest are significantly smaller than the horizon so that a 
Newtonian treatment is adequate. 

• Each element in the fluid expands to a maximum or turnaround radius (as measured from 
the center of the protogalaxy) before collapsing with inner regions reaching turnaround first. 
It is during the early stages of collapse that an outward moving shock develops. As infalling 
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material passes through the shock, it is heated rapidly and decelerated. 

• Vorticity is generated at the shock provided the velocity of the infalling gas is not every- 
where perpendicular to the shock surface. We demonstrate this for an axisymmetric prolate 
protogalaxy where the vorticity generated at the shock is in the azimuthal direction. 

• An external tidal torque applied to the protogalaxy generates shear angular momentum which 
in turn couples to the vorticity in the postshock region. As an example we consider again the 
model protogalaxy described above but now under the influence of a tidal torque along one 
of its short axes. The resultant shear field couples to the vorticity generated at the shock to 
yield a large-scale dipole-like vorticity field oriented along the direction of the tidal torque, 
i.e., along what will ultimately be the spin axis of the galaxy. The concomitant magnetic field 
has the same geometry and provides the seed field for subsequent dynamo action. 

Most previous analyses of protogalactic field generation have sought order of magnitude es- 
timates for the seed field strength without making direct contact to specific models of structure 
formation (Pudritz & Silk 1989; Lesch & Chiba 1994). Our results are in agreement with these es- 
timates and go one step further by providing a clear and simple picture of the geometry of the seed 
field. Our model is in the spirit of the semi-analytic and numerical studies of disk galaxy formation 
by Mestel 1963; Fall & Efstathiou 1980; Katz & Gunn 1991; Dalcanton, Spergel, & Summers 1997 
and others. However, in those works, angular momentum and vorticity are assumed ab initio and 
they are therefore unable to shed light on the creation of the first magnetic fields. In contrast, our 
model explicitly follows vorticity generation during the earliest stages of galaxy formation. 

In Section 2 we review the vorticity-magnetic field connection, derive an expression for the 
magnetic field generated at an oblique shock, and apply the results to our model protogalaxy. 
These analytic calculations are enough to obtain an estimate for the magnitude of the magnetic 
field as well as the general features of its geometry. The numerical simulations presented in Section 
3 provide a check of these results and also serve to illustrate some of the pitfalls inherent in using 
simulations to study problems of this type. These simulations do not include angular momentum 
and so in Section 4, we present a simple semi-analytic calculation for the postshock evolution of the 
vorticity and magnetic field in the presence of shear angular momentum. We conclude, in Section 
5, with a summary and discuss directions for future work. 

2. Vorticity Generation in a Protogalaxy: Analytic Treatment 

2.1. Vorticity-Magnetic Field Connection 

The evolution of a collisional fluid is described by the Euler equation 

3v 1 

— + (v-V)v = --Vp-VV> (1) 
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together with Poisson's equation for the gravitational potential ip, the continuity equation, an 
equation of state, and an energy equation. Taking the curl of eq. 1 yields the following for the 
vorticity u = V x v: 

du> „ , Vp x Vp 

— -Vx(vxc) = — -g 

Thus, while galaxies can acquire angular momentum through tidal fields, vorticity arises through 
purely gasdynamical processes, namely pressure and density gradients that are not colinear. 



V x (v x u) = (2) 



Biermann (1950) realized that a similar situation exists for magnetic fields. In the usual 
formulation of magnetohydro dynamics (MHD), the evolution of a magnetic field is described by 
the equation 

f -Vx(vxB)-A V =B = (3) 

where a is the conductivity. If B is initially zero, then it will be zero at all times. However, the 
derivation of eq. 3 assumes a form for Ohm's law that is not strictly valid for an electron-ion fluid. 
A careful treatment dictates that we include, on the right hand side, the term 

^ c Vn e x Vp e , 
r = Fe 4 

e nj 

where n e and p e are the number density and pressure of free electrons. Approximate local charge 

neutrality implies that n e ~ n p = xp/ m p where n p is the proton number density, \ ls the ionization 

fraction, and m p is the proton mass. In addition, since the electron temperature is expected to 

be approximately equal to the total gas temperature, p e ~ pn e / (n e + n p ) = px/ (1 + x). We can 

therefore write „ _ 

Vp x Vp 

T = a 5 (5) 

where a = m p c/e (1 + \) l-05x ~ 10~ 4 Gauss • s (Kulsrud et al. 1997). 

In the limit of vanishing diffusion, the equations for u> and B take identical forms. Together 

with the assumption that initially both the vorticity and magnetic field are zero, we have the 
relation, 

B = auj ~ 1(T 4 cj (6) 

where the units of B and u are Gauss and Hz respectively. The growth and evolution of the 
magnetic field therefore mirrors that of the vorticity up until the time when the diffusive effects of 
viscosity and conductivity become important 2 . 



2.2. Shock Wave Preliminaries 

For a barytropic fluid, p = p(p) and therefore T = 0. However, at curved shocks, the equation 
of state is more complicated (p = p(p, s) where s is the entropy) reflecting the fact that bulk kinetic 
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energy can be converted to thermal energy. We therefore expect r / and hence the generation 
of both vorticity and magnetic fields (Kulsrud et al. 1997). 

An ideal shock can be treated as a surface of discontinuity in the gas flow. By imposing certain 
jump conditions at this surface we can derive a relationship between the velocity field of the gas in 
the pre- and post-shock regions and hence an expression for the vorticity generated at the shock. 
In this way, we bypass eq. 5 and avoid dealing with the complicated gasdynamics that occurs inside 
the shock 

The standard shock wave jump conditions (e.g., Landau and Lifshitz 1997) consist of a set of 
three relations that guarantee the conservation of mass, momentum, and energy across the shock. 
Consider a point on a shock surface with velocity v s and normal n. Let po, vo, po, and ho be 
the density, velocity, pressure, and enthalpy in the preshock region and p, v, p, and h be the 
corresponding quantities in the postshock region. The jump conditions are 

pu = p u (7) 
p + pu 2 = p + p ul (8) 

h + ^u 2 = h + -u 2 (9) 
where no and u are the velocity components along n in the rest frame of the shock: 

= (vo - v s ) • n u = (v - v s ) • n . (10) 

In addition, we have that the component of the velocity tangent to the shock is continuous. 

For the situation at hand, the preshock gas is relatively cold and we can therefore set po ~ 
~ /to- The three jump conditions are then easily combined to give 

where 7 is the polytropic index of the infalling gas (equal to 5/3 for an ideal gas). This allows us 
to express the velocity in the postshock region in terms of the preshock velocity: 

v = v + (/( 7 ) - l)n n (12) 

where 7(7) = (7 - 1)/ (7 + 1). 



2.3. Vorticity in the Postshock Region 

The velocity field in the postshock region is determined not only from the initial velocity field 
and the geometry of the shock (through eq. 12) but also from the evolution of the gas once it has 
passed through the shock. All of this is incorporated into numerical simulations discussed in the 
next section. Here we present an analytic model for an idealized protogalaxy. The key simplification 
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is to ignore the evolution of the gas in the postshock region. Our picture is that the shock wave 
sweeps through the gas transforming the velocity field from vo to v according to eq. 12 where uo 
and n are evaluated at each point at the instant when the shock passes through. In this way, no 
and n can be treated as functions of position and we can then calculate the vorticity by taking the 
curl of v. As discussed above, we expect v to be curl-free so that 

w = (/( 7 )-l)Vx( Uo fi) (13) 

To proceed further we require a specific model protogalaxy. As a starting point, we consider 
the spherical infall model (e.g., Gunn 1975, Gott 1977, Fillmore & Goldreich 1984, Bertschinger 
1985, Ryden & Gunn 1987, Ryden 1988) wherein matter is divided into spherical shells which 
expand to a maximum or turnaround radius and then collapse toward the center. In the case of 
collisional matter (e.g., Bertschinger 1985) infalling shells are decelerated and heated as they pass 
through an outward moving shock. 

Current theories of structure formation present a far more complicated picture than that rep- 
resented by the spherical infall model. In particular, structure formation is believed to proceed 
hierarchically, with subgalactic objects forming first, and then coalescing to form galaxies and clus- 
ters. For our purposes, the key deficiency of the spherical infall model is its restriction to spherical 
symmetry, since this precludes vorticity generation. Moreover, a spherically symmetric protogalaxy 
cannot acquire angular momentum through tidal torques. Of course, there is no reason to expect 
protogalaxies to be spherically symmetric. Indeed, the halos found in collisionless N-body simu- 
lations are generally triaxial, with prolate shapes favored slightly over oblate ones. Furthermore, 
the angular momentum vector for these systems is generally aligned with the short axis of the halo 
(Carlberg & Dubinski 1991; Warren et al. 1992). These results motivate us to consider a simple 
model protogalaxy that forms from an axisymmetric, prolate density perturbation. An external 
tidal torque, applied perpendicular to the symmetry axis of the perturbation, generates shear an- 
gular momentum. In the spirit of the spherical infall model, we assume that the perturbation is 
smooth and featureless with a density profile that decreases with radius. The protogalaxy therefore 
forms from the inside out. 

For the moment, we ignore tidal torques. The perturbation will therefore evolve into an 
axisymmetric protogalaxy. Consider a spherical coordinate system, (r, 0, (ft), with polar axis oriented 
along the symmetry axis of the protogalaxy. The density and pressure gradients that occur, for 
example, at an outward moving shock, will be in the f and 9 directions which imply that any 
vorticity generated will be along the cj) direction. Moreover, by symmetry the vorticity above and 
below the equatorial plane will be in opposite directions. This is not surprising since, in the absence 
of tidal fields, the net circulation of the system must be zero. 

It is instructive to consider a simple ansatz for the evolving protogalaxy. Specifically, we 
assume that isodensity contours for the gas are concentric spheroids, i.e., p = p(r s ) where f 2 = 
r 2 (sin 2 9 + cos 2 9/q 2 ). q is the flattening parameter which, in general depends on time and radius. 
For simplicity, we will ignore this complication (we focus on a small region in the neighborhood of 
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the shock) and further assume that the deviation from spherical symmetry is small (\q — 1| <C 1). 
The shock surface is described by the equation 

f2(t) = r 2 (sin 2 # + cos 2 fl/g 2 ) (14) 

and the normal to this surface is, to first order in (q — 1), given by 

n = r + 2(g-l)sin<9cos6>0 (15) 

Thus, one contribution to uj will be of the form 2 (/(7) — 1) (1 — 5) (uo/f s ) sin# cos 00. There is a 
second contribution that is proportional to Viio x n which has a similar form 3 and we can therefore 
write 

u oc (q — 1) (uo/r s ) sin cos 6*0 (16) 

The magnitude of the vorticity is therefore set by the velocity of the infalling gas (in the rest frame 
of the shock) divided by the shock radius. This is roughly equal to the reciprocal of the turnaround 
time, T ta , a result we might have anticipated from dimensional analysis. As a specific example, 
consider Bertschinger's (1985) solution for secondary infall of collisional matter onto an already 
collapsed overdensity. In this self-similar model, the turnaround radius at time t is r ta (t) oc t 8 / 9 
and the radius of the shock is r s (t) = A s r ta (i) where A s is a constant ~ 0.33 for 7 = 5/3. Likewise 
the velocity of the gas, immediately before passing through the shock, is given by vq = —Vr ta /t 
where V ~ 1.47 for 7 = 5/3. This implies that u /r s ~ 1.7/T ta . 

For a galaxy-sized object, T ta ~ 10 16 s and therefore uj^ ~ 10~ 16 s _1 . This is roughly a factor of 
10 less than the local value of the vorticity in the Milky Way as determined from the Oort constants 
(e.g., Binney & Merrifield 1998), a reasonable result given that the vorticity will be amplified during 
the formation of the disk itself. The strength of the corresponding magnetic field is ~ 10~ 20 Gauss. 
We will return to this result in the next section. 



3. Numerical Simulations 

In this section we present the results of numerical simulations that are designed to test and 
augment the analytic model described above. The simulations follow the evolution of an isolated 
axisymmetric density perturbation in an otherwise flat (Einstein-de Sitter) universe. The cos- 
mic fluid consists of dark matter and gas in a 10:1 ratio. The simulations are performed using 
HYDRA (Couchman, Thomas, &; Pearce 1995): Gravitational forces are calculated with an adap- 
tive particle-particle particle-mesh (AP 3 M) algorithm while gasdynamics is treated using smooth 
particle hydrodynamics (SPH). Simulations are run with 32 3 particles of each species. 



3 To evaluate this term, we require an ansatz for the preshock gas flow, vo. Since this is assumed to be irrotational, 
it can be written as the gradient of a scalar function. A reasonable ansatz (akin to the Zel'dovich approximation) is 
vo oc Vi/>. In any case, we expect that vo -0/vo ■? = 0(\q— 1|) and likewise for v s so that uo = vo r —v ar + 0((q~ l) 2 ). 
Moreover, uo should have the form uo(r, 0) — Uq + u%{q — 1) cos 2 (9) where uj and Uo are functions of r which depend 
on the details of the model. 
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The initial density profile has the form p(r, 0) = Pfe(i) (1 + o{f)) where pb(t) is the background 
density for an Einstein-de Sitter Universe, r = r (sin 2 + cos 2 0/q 2 ) 1 ^ 2 , and 



<5(r) 



^o(l-^2^) 2 ) r<r, 

i* (* fen 



(17) 

r r <r < R 



To set up initial conditions, we begin with an interlaced lattice of gas and dark matter particles. 
Those particles a distance R from a chosen center are discarded and the ones that remain arc 
displaced from their original lattice sites so as to achieve the desired density profile. Velocities are 
then assigned according to the Zel'dovich approximation. In the simulations presented here, 5q = 1, 
r c /R = 0.4, a = 2, and for the prolate runs, q = 1.5. With this choice of parameters, the mean 
density enhancement is 5 = 0.377. 

The simulation units are such that the total mass M = 1 and Newton's constant G = 0.0194. 
Neither cooling nor star formation are included in the simulations and so there is some freedom 
in choosing units for dimensional quantities. In conventional models of structure formation, such 
as the Cold Dark Matter scenario (CDM) and its variants, 6 can be identified with the rms mass 
fluctuation on a scale R: 

a M (t) ^ ((AM/M) 2 ) 1 / 2 = (i + ^ —- P{k)W \kR)j (18) 

where M ~ 1.2 X lO 12 /i 2 M (_R/Mpc) 3 is the total mass in a sphere of radius R, P(k) is the 
linear power spectrum for the model, h is the present value of the Hubble parameter in units 
of lOOkms^Mpc- 1 , z is the redshift and W(x) = 3(sinx — xcosx) /x 3 is the top hat window 
function. Thus, by setting gm = 0.377, we can determine the initial redshift for the simulation, 
Zi, as a function of M. It is then straightforward to relate simulation units to physical units. This 
is done, for three representative masses, in Table 1 where, in computing gm, we have assumed a 
spatially flat CDM universe with h = 0.7, 0,b = 0.05, and COBE normalization (cr 8 ~ 1.7). The 
transfer function was calculated using the fitting formula of Eisenstein and Hu (1999). 

Figure 1 presents the phase space particle distribution (radius r vs. radial velocity v r ) for the 
spherical run. We see that the turnaround radius, r ta , increases with time. For r > r ta , the gas 
and dark matter particles evolve as a single fluid. For r < r ta , the dark matter particles exhibit 
multiple phase space streams that are characteristic of collisionless infall (cf Figure 10 of Fillmore & 
Goldreich 1984 and Figure 6 of Bertschinger 1985). Note however that at late times, these streams, 
especially the outward moving ones, become rather chaotic. This is a result of an instability in the 
spherical infall model first described by Henriksen k, Widrow (1997). In contrast, the gas particles 
are decelerated with v r — > for r — > 0. The presence of an outward moving shock is clearly seen in 
Figure 2 where we plot the temperature T as a function of r for frames (b) and (d) of Figure 1. 

We next turn to vorticity. In SPH one determines the evolution of a fluid system by following 
the motion of fiducial particles which are labeled with local kinematic and thermodynamic quantities 
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(for a review, see Monaghan 1992 and references therein). Any function /(r) of these quantities 
may be approximated by the following summation: 

N 

f{ Y ) = Y J —W{T-n;R)f b (19) 
where /& is the value of /(r) for the 6'th particle and 

N 

Pb = J2^bW(r-r b ;R) (20) 
6=1 

In these expressions, W is a user-supplied window function with characteristic radius R. This 
measurement process introduces an error which can be minimized (but not eliminated) by an 
appropriate choice of R. Quantities that involve gradients require a bit more care and the measure- 
ment prescription is not always unique. Following Monaghan (1992) the vorticity is determined as 
follows: 

1 N 

u a = (V x v) a = — m bv ba x V a W ab (21) 

Pa 6=1 

where v& a = v& — v a , W ab = W(r a — r^,; R) and ^afflab denotes the gradient of W ab with respect 
to r a . Following the usual practice, we choose the r a to be the positions of the particles themselves 
though in principle r a can be taken to be at any point in the simulation volume. 

In addition to the measurement error discussed above, there is an error associated with the 
integration of particle orbits. Moreover, the "boxy" nature of the particle distribution, an artifact 
of the initial conditions, will lead to vorticity generation even with q = 1 since we do not have true 
spherical symmetry. 

As a diagnostic test of these potential difficulties we determine the vorticity field in the spherical 
run discussed above. The result is shown in Figure 3 where we plot the (r, 9, <p) components of the 
measured a; as a function of r for frame (b) of Figure 1. The large 6 and <p components imply 
that there are angular gradients in v r and/or radial gradients in vq and v<f>. For exact spherical 
symmetry and properly treated gasdynamics, the velocity fields should be purely radial and the 
only gradients in the r direction, ujq and uj^ represent the first terms that arise when spherical 
symmetry is broken. In contrast, a nonzero u r requires angular gradients in the tangential velocity 
field and is therefore second order in small quantities and so it is not surprising that the measured 
uj r is the smallest of the three components. As expected, the amplitudes of ujq and oj^ decrease 
with increasing particle number, roughly as iV~ 1//3 . 

As mentioned above, errors are introduced into the vorticity calculation simply because we 
are attempting to determine a continuous field from information at discrete and irregularly spaced 
points. In order to quantify this aspect of the problem, we calculate the vorticity for a distribution 
of particles with the same positions as those used to generate Figure 3 but with velocities cho- 
sen by hand to reproduce a prescribed velocity field. Equation 21 is then used to determine the 



-10- 



"measured" vorticity field. For this experiment, we assume a prescribed velocity field of the form 

1 /2 

v z = (x 2 + y 2 ) which implies a constant vorticity field, oj = 4>. The measured field, shown in 
Figure 4, indicates that for most particles, the SPH prescription does a good job of calculating the 
vorticity. However, for a subset of particles, errors of order unity are introduced. 

The difficulties inherent in following the generation and evolution of vorticity in hydrodynamic 
simulations is apparent in the simulations of Kulsrud et al. (1997). They determine the magnetic 
field by solving eq. 3 with the additional term eq. 5 included on the right hand side. As a check, they 
compare the result with the vorticity (scaled by the appropriate constant) and find discrepancies 
of order unity (cf. their Figure 4). 

We next calculated the vorticity for the prolate protogalaxy simulation (Figure 5). As expected, 
there is now significant vorticity generated at the shock, primarily in the azimuthal direction. The 
vorticity in the r and 9 directions is again an artifact of the simulation. The rms of lo^ is a factor 
of 5 greater than that of ujq and uj t . This can be viewed, in some sense, as a measure of the 
signal-to-noise of the simulation. 

In Figure 6, uj^ is plotted as a function of 9. The expected antisymmetry about the equatorial 
plane {9 = vr/2) is readily apparent. In particular, the vorticity near the poles (9 = and ir) 
vanishes. However, we also find that there is vorticity generated with the "wrong sign", i.e., uj^ < 
for < 9 < 7r/2 and lo^ > for tt/2 < 9 < tt. This can be understood as follows: When the gas 
flows through the shock, it is refracted away from the symmetry axis and toward the equatorial 
plane. This leads to a region of high pressure and density in the equatorial plane forcing the gas to 
move out along the symmetry axis and creating a region of vorticity with the opposite sign. This is 
illustrated in Figures 7a and 7b where we show the velocity field of the particles in the simulation. 

For a system mass of 7 x 1O 11 M , corresponding to a spiral galaxy roughly the size of the Milky 
Way. The magnitude of the vorticity is ~ 10~ 15 s _1 in good agreement with our earlier estimate. 
The magnitude of the corresponding magnetic field is ~ 10~ 19 G. The protogalaxy at these early 
stages is roughly 25 kpc in size whereas the actual disk will have a radius ~ 10 kpc and a thickness 
~ 1 kpc. Contraction of the protogalaxy in the plane of the disk will therefore amplify the seed field 
by a factor of (25/10) 2 ~ 6 while collapse perpendicular to this plane will amplify the field by a 
factor ~ 25 (Lesch & Chiba 1995). We therefore expect a field strength in the fully assembled disk 
galaxy of 1.5 x 10 -17 Gauss. This is approximately what is required to seed the galactic dynamo. 

An alternative scenario is to generate the first magnetic fields in 10 6 M objects. The seed 
fields are approximately two orders of magnitude larger (see Table 1). More importantly, the 
dynamical time for these systems is significantly shorter. It may therefore be possible for dynamo 
action to amplify fields on these scales before the disk is assembled. 
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4. Post-Shock Evolution 

In the axisymmetric model described above, the vorticity and magnetic field generated at the 
shock are in the azimuthal direction and are antisymmetric about the equatorial plane. Mixing 
of gas from above and below this plane will lead to a rapid decrease in the vorticity, a reflection 
of the fact that angular momentum has not been included. The evolution of the magnetic field 
involves recombination and is therefore more complicated. It is however clear that no large-scale 
coherent field will survive without the addition of angular momentum. As discussed above, shear 
angular momentum is generated by tidal interactions with neighboring protogalaxies and is typically 
oriented along one of the short axes of the protogalaxy. It is the action of the shear field on the 
vorticity and magnetic fields that leads ultimately to a dipole configuration for these fields. This 
process can be illustrated by the following simple calculation. A set of particles are used to represent 
fluid elements labeled by their position, velocity, velocity gradient, and magnetic field. We assume 
force-free evolution so that each particle evolves independently according to the following (Cartesian 
coordinate) equations: 

dxi dvi 

— — = Vi — — = 22) 

dt dt v ; 

= (diV k ){d kVj ) ^ = BjdjVi- Bid jVj (23) 

Initially, the magnetic field is in the azimuthal direction and is antisymmetric about the equatorial 
(xz) plane (Figure 8a). The prescribed velocity field includes shear angular momentum about the 
z-axis and an inward radial flow. The latter is meant to model the continual contraction of the 
fluid under the influence of gravity. After a short period of time, these field lines are sheared into 
a dipole configuration (Figure 8b). 



5. Conclusions 

In summary, we have presented a detailed investigation of magnetic field generation during the 
collapse of a protogalactic density perturbation. The first fields appear, via the Biermann battery 
effect, in the region of the outward moving shock that develops in the collapsing protogalaxy. Shear 
angular momentum is then able to reconfigure the field into a dipole pattern oriented along the 
spin axis of the protogalaxy. The predicted field strength, once the disk has formed, is estimated 
to be 10~ 17 Gauss. With a seed field of this magnitude, dynamo action can create microgauss fields 
by the current epoch. The magnitude of the associated vorticity at the time of disk formation is 
roughly equal to its present day value. 

Virtually all galactic dynamo models assume azimuthal symmetry with respect to the spin axis 
of the disk. The dynamo equations in these models possess an invariance with respect to reflections 
about the equatorial plane and therefore the solutions can be divided into two groups: odd modes 
which consist of a dipole-like poloidal field together with an antisymmetric toroidal field and even 
modes which consist of quadrupole-like poloidal fields together with symmetric toroidal fields. In 
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general, even modes are favored (faster growing) when the fields are confined to a disk, while odd 
modes are favored in more spherical configurations (e.g., Ruzmaikin, Shukurov, & Sokoloff 1988). 
Observations would seem to indicate that both types of configurations are present in nature: The 
fields in the inner regions of the Milky Way, for example, appear to be predominantly antisymmetric 
with respect to the disk plane (Han et al. 1997) while those in M31 are evidently symmetric (Han, 
Beck, & Berkhuijsen 1998). 

Our analysis suggests that dipole-like seed fields are favored (see, also Krause & Beck 1998). 
However, in a more realistic model, based on hierarchical clustering, we expect both dipolar and 
quadupolar fields to be produced. 

If galactic magnetic fields have their origin in the Biermann battery effect operating in pro- 
togalactic shocks, it should be possible to follow the formation of a disk galaxy from primordial 
density perturbation, with B = 0, to mature galaxy with a microgauss field. Our analysis and 
simulations have taken the first step in this ambitious program. Future work will include cosmo- 
logical tidal fields as well as small scale perturbations. In addition, the magnetic field will have to 
be treated explicitly since the correspondence with vorticity is ultimately lost. 

We wish to thank R. Henriksen, H. Couchman, and, L. Chamandy for useful discussions. LMW 
acknowledges the hospitality of The University of Chicago during a sabbatical stay and GD ac- 
knowledges the hospitality of the Canadian Institute for Theoretical Astrophysics. This work was 
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Canada. 



-13- 



REFERENCES 

Bertschinger, E. 1985, ApJS, 58, 39 
Biermann, L. 1950, Zs. Naturforsch., 5a, 65 

Binney, J. J. and Merrifield, M. 1998 Galactic Astronomy (Princeton, NJ: Princeton University 
Press) 

Bisnovatyi-Kogan, G. S., Ruzmaikin, A. A., and Syunyaev, R. A. 1973, Sov. Ast. 17, 137 

Carlberg, R. G. and Dubinski, J. 1991, ApJ, 369, 13 

Couchman, H.M.P., Thomas, P. A. , Pearce, F. R. 1995, ApJ, 452, 797 

Dalcanton, J. Spergel, D., and Summers, F. 1997, ApJ, 482, 659 

Eisenstein, D. J. and Hu, W. 1999, ApJ, 511, 5 

Fall, S. M. and Efstathiou, G. 1980, MNRAS, 193, 189 

Field, G. B. and Carroll, S. M, 1998, astro-ph/9811206 

Fillmore, J. A. and Goldreich, P. 1984, ApJ, 281, 1 

Gott, J.R. 1975, ApJ, 201, 296 

Gunn, J. E. 1977, ApJ, 218, 592 

Han, J. L., Manchester, R. N., Berkhuijsen, E. M., and Beck, R. 1997, A&A, 322, 98 
Han, J. L., Beck, R., and Berkhuijsen, E. M., 1998, A&A, 335, 1117 
Harrison, E. R. 1970, MNRAS, 147, 279 
Harrison, E. R. 1973, Phys. Rev. Lett., 30, 18 

Henriksen, R.N. and Widrow, L.M. 1997, Phys. Rev. Lett., 78, 3426 

Hoyle, F. 1949, in Proc. Symposium on Motion of Gaseous Masses of Cosmical Dimensions, Prob- 
lems of Cosmical Aerodynamics, ed. J. M. Burgers and H. C. van de Hulst (Dayton: Central 
Air Documents Office), p. 195 

Katz, N. and Gunn, J. E. 1991, ApJ, 377, 365 

Kim, E.-J., Olinto, A. V., and Rosner, R. 1996, ApJ, 468, 28 

Krause, F. and Beck, R. 1998, A&A, 335, 789 

Kronberg, P. P. 1994, Rep. Prog. Phys. 325 



-14- 



Kulsrud, R.M., Cen, R., Ostriker, J.R and Ryu, D. 1997, 480, 481 

Landau, L. D. and Lifshitz, E. M. 1959, Fluid Mechanics (London: Pergamon Press) 

Lesch, H. and Chiba, M. 1995, A&A , 297, 305 

Mestel, L. 1963, MNRAS, 126, 553 

Monaghan, J. J. 1992, ARA&A, 30, 543 

Peebles, P. J. E. 1969, ApJ, 155, 393 

Pudritz, R.E. and Silk, J. 1989, ApJ, 342, 650 

Quashnock, J., Loeb, A. and Spergel, D. 1989, ApJ, 344, L49 

Ratra, B. 1992, ApJ, 391, LI 

Rees, M. J. QJRAS, 1987, 28, 197 

Ryden, B. S. 1988, ApJ, 329, 589 

Ryden, B.S. and Gunn, J.E. 1987, ApJ, 318, 15 

Ruzmaiken, A. A., Shukurov, A. M., and Sokoloff, D. D. 1988, in Astrophysics and Space Science 
Library, Magnetic Fields in Galaxies (Dordrecht: Kluwer) 

Sicotte, H. 1997, MNRAS, 287, 1 

Turner, M. S. and Widrow, L. M. 1988, Phys. Rev. D37, 2743 

Warren, M. S., Zurek, W. H., Quinn, P. J., and Salmon, J. K. 1992, ApJ, 399, 405 
Wasserman, I. 1978, ApJ, 224, 337 
White, S. D. M. 1984, ApJ, 286, 38 
Vachaspati, T. 1991, Phys. Lett. B, 265, 258 

Zel'dovich, Ya. B., Ruzmaikin, A. A., and Sokoloff, D.D. 1983, Magnetic Fields in Astrophysics, 
Gordon and Breach, New York 



This preprint was prepared with the AAS IATgjX macros v5.0. 



-15- 



Table 1. Physical units for numerical simulations 



[M] Mq 


10 6 


10 9 


7 x 10 


Zi a 


120 


63 


25 




43 


22 


8.4 


[L] (kpc) 


5.9 


11 


24 


[V] (kms- 1 ) 


14 


100 


560 


u] (lO-^s- 1 ) 


75 


29 


7.5 


[B] (10- 20 G) 


79 


31 


7.9 



a Redshift at the start of the simulation 
a Redshift corresponding to panel (d) of 
Figure 1 
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Fig. 1. — Phase space (radius r vs. radial velocity v r ) distribution of particles in the spherical run. 
The blue points represent dark matter particles while the red points represent gas particles. 
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Fig. 2. — Temperature as a function of radius corresponding to frames (b) and (d) of figure 1 
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Fig. 3. — (r,6,<p) components of u for frame (b) of Figure 1. 
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Fig. 4. — Same as Figure 3 with a prescribed velocity field corresponding to u> = (j>. 
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Fig. 5. — Same as Figure 3 but for the prolate protogalaxy (q = 1.5). 
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Fig. 6. — W0 as a function of 0. 
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Fig. 7. — Velocity field for a prolate protogalaxy (q = 1.25) as a function of cylindrical radius r and 
position along the symmetry axis z. (a) The entire protogalaxy. (b) Inner region of the protogalaxy 
corresponding to the rectangular box in (a). 
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Fig. 8. — Magnetic field in the Lagrangian calculation described in Section 4. (a) The initial 
field configuration. The symmetry axis of the protogalaxy coincides with the y-axis. (b) Final 
configuration. Particles have evolved under the influence of a shear field that has net angular 
momentum in the z direction. 



